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Fermi liquid theory, the standard theory of metals, has been challenged by a 
number of observations of anomalous metallic behavior found in the vicinity 
of a quantum phase transition. The breakdown of the Fermi liquid is accom- 
plished by fine-tuning the material to a quantum critical point using a control 
parameter such as the magnetic field, pressure, or chemical composition. Our 
high precision magnetization measurements of the ultrapure /-electron based 
superconductor /3-YbAlB 4 demonstrate a scaling of its free energy indicative 
of zero-field quantum criticality without tuning in a metal. The breakdown of 
Fermi-liquid behavior takes place in a mixed- valence state, in sharp contrast 
with other known examples of quantum critical /-electron systems that are 
magnetic Kondo lattice systems with integral valence. 
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Quantum phase transitions occur at zero temperature as a consequence of quantum, rather 
than thermal correlations. Generally, a quantum critical point (QCP) can be reached by driving 
a finite temperature critical point (1, 2) or a first-order critical end-point (3) to absolute zero. The 
breakdown of Fermi liquid (FL) behavior in metals observed near a magnetic QCP challenges 
our current understanding of strongly correlated electrons. While the mechanism of unconven- 
tional quantum criticality is actively debated, there is a growing consensus that the underlying 
physics involves a jump in the Fermi surface volume associated with a partial electron localiza- 
tion (4-8). To date, the FL breakdown has only been observed by fine-tuning a material to a 
QCP using a control parameter such as magnetic field, pressure, or chemical composition. 

Recent work, reporting the discovery of superconductivity in an ytterbium-based heavy 
fermion material /3-YbAlB 4 , raised the interesting possibility that this system may exhibit quan- 
tum criticality without tuning (9). In this compound, signatures of quantum criticality were 
observed to develop above a tiny superconducting (SC) dome, with a SC transition temperature 
of T c ~ 80 mK and an upper critical field [IqB.^ ~ 30 mT (9, 10). While this observation moti- 
vated the possibility of a zero field quantum critical point, it did not rule out a QCP located near 
the upper critical field H c2 , as observed in the heavy fermion superconductor CeCoIn 5 (11). 

In this report, we present clear and quantitative evidence that quantum criticality develops at 
zero field without tuning in /3-YbAlB 4 , buried deep inside the SC dome. Moreover, we report 
a simple T/5-scaling form of the free energy spanning almost four decades in magnetic field, 
revealing that the signatures of the putative quantum critical point extend up to temperatures T 
and fields B more than 100 times larger than T c and HqH C 2, respectively. 

To quantify the free energy F(T,B), we employed high-precision measurements of the 
magnetization M = —dF/dB. Measurements were made on ultra high purity single crystals 
with a mean free path exceeding 1000 A and residual resistivity less than 0.6 iiVt cm, which 
were carefully etched to fully remove surface impurities (12). Our measurements revealed a 
simple T / B scaling over a wide range of temperature and field, governed by single quantum- 
critical (QC) scaling exponent previously masked (9) by a limited experimental resolution and 
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the impurity effects caused by surface and bulk impurities (72). The T/B scaling leads to the 
following significant consequences. First, the quantum critical physics is self-similar over four 
decades of T/B, with no intrinsic energy scale. Second, the field-induced Fermi liquid is char- 
acterized by a Fermi temperature that grows linearly with the field, determined by the Zeeman 
energy of the underlying critical modes. Finally, the scaling allowed us to determine an upper 
bound on the magnitude of the critical field \B C \ < 0.2 mT, which is well inside the SC dome 
and comparable with the Earth's magnetic field: this indicates that /3-YbAlB 4 is intrinsically 
quantum critical, without tuning the magnetic field, pressure, or composition. 

These results are surprising, given the fluctuating valence nature of this material, with va- 
lence Yb +2 75 significantly far from integral, revealed by recent experiments (13). All QC 
heavy-fermion intermetallics known to date have an almost integral valence which stabilizes 
the local moments (1, 2). Such so-called Kondo lattice systems are characterized by a small 
characteristic scale T , below which the moments are screened to form a paramagnetic heavy 
FL. Various types of order, such as superconductivity and antiferromagnetism (AFM), com- 
pete with the heavy FL, leading to quantum criticality, as seen in, for example, CeCu5.gAu .i 
(T = 6.2 K) (14) and YbRh 2 Si 2 (24 K) (2). By contrast, mixed-valence compounds display 
a much larger T , below which they typically behave as stable FLs with moderate quasiparti- 
cle effective masses and no competing order. For example, YbAl 3 , with non-integral valence 
Yb +2 n , is characterized by T ~ 300 K (15). 

A remarkable feature of /3-YbAlB 4 (Fig. 1A) is that it is quantum critical (9), yet the scale 
T ~ 250 K, obtained from the resistivity coherence peak, is one or two orders of magnitude 
larger than in other known QC materials. This is confirmed by the scaling behavior of the mag- 
netic specific heat: Cm/T = |? In (T /T), where So is a constant (inset of Fig. IB). The — In T 
dependence of Cm/T in the three QC materials CeCus.gAuo.i (14), YbRh 2 Si2 (2, 16), and (3- 
YbAlB 4 collapse onto one curve after setting T for /3-YbAlB 4 ~ 200 K. The recent observation 
of intermediate valence (Yb +2 75 ) in /3-YbAlB 4 at 20 K using hard X-ray photoemission spec- 
troscopy (13) is consistent with this large T . 
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The quantum criticality in this valence fluctuating state is accompanied by several properties 
reminiscent of an integral- valence Kondo lattice. To understand their origin, it is useful to com- 
pare them with those of «-YbAlB 4 (Fig. 1A) (17), which is a locally isostructural polymorph 
of /3-YbAlB 4 and a FL with a similarly intermediate valence (Yb +2 73 ) (13). Instead of Pauli 
paramagnetism normally seen in a valence fluctuating material, the magnetic susceptibilities of 
both materials display Curie- Weiss behavior with Weiss temperature 6w, X = C/(T + 8w), 
indicating the existence of local moments (Fig. 2A). In addition, both materials have a maxi- 
mum in —dM/dT at T* ~ 8 K, signaling a crossover from local moment behavior (Fig. SI). 
Below T*, Cm/T of the a phase levels off to a constant characteristic of heavy FL behavior, 
whereas that of the (3 phase continues to diverge (Fig. IB). Thus, the fate of local moments 
found above T* is different in these locally isostructural systems: Yb spins are fully screened 
in the a phase, but may well survive down to lower temperatures in the phase and produce 
the quantum criticality. In both phases, strong correlation effects are manifest, for example, in 
the strongly enhanced ^| T ^ ~ 130 mJ/mol K 2 , two orders magnitude larger than the band 
calculation estimate (~ 6 mJ/molK 2 ) (18). 

These signatures indicate that both phases are governed by two distinct energy scales: a 
high-energy valence fluctuation scale T ~ 200 K, and a low-energy scale T* ~ 8 K, character- 
izing the emergence of Kondo-lattice physics. A possible origin of this behavior is the presence 
of ferromagnetic (FM) interactions between Yb moments, manifested by the large Wilson ratio 
i?w between \ an d Cm/T, observed in both a- and /3-phases (9) (minimum estimate R w > 7), 
and further corroborated by the observation of an ESR signal (19), generally only seen in the 
presence of FM correlations (20). Such FM interactions are known to give rise to Kondo reso- 
nance narrowing (21) in rf-electron systems where Hund's coupling causes a marked reduction 
in the Kondo temperature (21, 22). 

In the present case, the role of Hund's coupling is played by FM intersite RKKY interac- 
tions, probably along the short Yb-Yb bonds that form chains along the c-axis. The moments 
of a few n neighboring Yb ions may thus become aligned, forming a fluctuating "block" spin 
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S = nJ. The observed valence Yb +2 75 could then be understood in terms of Yb 3+ ^Yb 2+ 
fluctuations, as at any one time, approximately 1/4 of the Yb atoms along the chains are in a 
singlet Yb 2+ configuration, forming the ferromagnetic blocks of approximately n ~ 3 spins. 
The effect of these block spins is to exponentially suppress the characteristic spin fluctuation 
scale (21), resulting in localized-moment behavior. The absence of long-range magnetic order 
in the a or phase points to the presence of competing magnetic interactions. Indeed, the Weiss 
temperature 9\y ~ — 110 K (Fig. 2A), characteristic of an AFM, indicates the importance of 
magnetic frustration. The competing interplay of FM interactions and valence fluctuations thus 
leads to Kondo-lattice-like behavior in a mixed valent material, setting the stage for quantum 
criticality to emerge at lower temperatures in /3-YbAlB 4 . 

A prominent feature of the quantum criticality in /3-YbAlB 4 is the divergence of the mag- 
netic susceptibility x as T — > 0. By examining the field evolution of magnetization M = 
—dF/dB as a function of both T and B (12), we can accurately probe the free energy F near 
quantum criticality. Figure 2A shows the T dependence of x(B) = M / B for different values of 
B || c. Spanning four orders of magnitude in T and B, the data show a systematic evolution 
from a non-Fermi liquid (NFL) metal with divergent susceptibility at zero field (x ~ T" 1 / 2 ) to 
a FL with finite x ma field g^B > k B T. 

Intriguingly, the evolution of M/B found in the region T < 3 K and B < 2 T (see the inset 
of Fig. 2A) can be collapsed onto a single scaling function of the ratio T/B: 



as shown in Fig. 2B. The peak of the scaling curve lies at ksT / g\i?>B ~ 1, marking a cross-over 
between the FL and NFL regions, showing that ksTp ~ g^B plays the role of a field-induced 
Fermi energy, as shown in the inset of Fig. 2B. Integrating both parts of Eq. ([I]), one obtains the 
following scaling law for the free energy (12): 



where / is a scaling function of the ratio T/B with the limiting behavior: f(x) oc x 3 ^ 2 in the 
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(1) 




(2) 



NFL regime (x ^> 1) and f(x) oc const + x 2 in the FL phase (x <C 1). Indeed, the observed 
scaling of dM/dT in Eq. (1) is best fitted with <f>(x) = Ax(A + x 2 )^~ 2 (12), resulting in a 
particularly simple form of the free energy: 

*bc = -TT^h ((9^B) 2 + (k B T) 2 f 4 , (3) 
(k B T) l -/ z 

with the best fit obtained with effective moment g/j, B = 1.94/xb an d the energy scale k B T « 
6.6 eV of the order of the conduction electron bandwidth (12). This means that the free en- 
ergy depends only on the distance from the origin in the (T,B) phase diagram, similarly to the 
T/B scaling established in the Tomonaga-Luttinger liquids in one-dimensional metals (23, 24). 
Equation |3] implies that the effective mass of the quasi-particles diverges as m* ~ B~ l l 2 at 
the QCP (12). This divergence in a 3D material, together with the T/B scaling, cannot be ac- 
counted for by the standard theory based on spin-density-wave fluctuations (25, 26). Instead, it 
indicates a breakdown of the FL driven by unconventional quantum criticality. 

The T/B scaling suggests that the critical field B c of the quantum phase transition is actually 
zero. A finite B c would require that the argument of the scaling functions f(x) and <p(x) is the 
ratio x = T/\B — B c \, as seen for instance in YbRh 2 Si 2 (16). To place a bound on B c , we 
substituted this form for x into Eq.Q], seeking the value of B c that would best fit the experimental 
data. The Pearson's correlation coefficient R obtained for this fit (inset of Fig. 2B), indicates 
that B c is optimal at —0.1 ± 0.1 mT. The uncertainty is only a few times larger than the Earth's 
magnetic field (~ 0.05 mT). More significantly, it is two orders of magnitude smaller than 
/i i?c2 = 30 mT, and strikingly six orders of magnitude smaller than valence fluctuation scale 
T ~ 200 K. Thus /9-YbAlB 4 provides a unique example of essentially zero-field quantum 
criticality. 

Further evidence for zero-field quantum criticality is obtained from an analysis of the mag- 
netocaloric ratio, IV = — ^H^f = — 9M ^ 9T (pig 3) Here, C is the total specific heat (12). 
Our results show a clear divergence of Th/B as T — > in the NFL regime, which is a strong 
indicator of quantum criticality (27). From the NFL regime, we can extract the critical field 
\B C \ < 0.2 mT, consistent with the estimate of B c obtained from the scaling behavior of M, 



Eq.dRFig. S2). 

The remarkably simple T/B scaling in the thermodynamics enables us to characterize the 
QC excitations of /3-YbAlB 4 . In particular, the collapse of all magnetization data in terms of 
the dimensionless ratio r = k^T / (g^B) between the Boltzmann energy fcgT and the Zeeman 
energy g^B indicates an absence of scale in the zero-field normal state. Furthermore, the 
appearance of a field-induced Fermi energy, linear over more than three decades in B, shows 
that the underlying critical modes are magnetic in character. 

Using the Heisenberg energy-time uncertainty principle (AtAE>h), we can reinterpret the 
T/B scaling in the time domain, visualizing the field-induced FL as a kind of "quantum soda" of 
bubbles of quantum critical matter of finite duration tq = h/ g^B, immersed in a FL. At finite 
temperatures, thermodynamics averages the physics over a thermal time scale tt = h/k B T, 
thus the quantity r = tq/t t ~ T/B in the scaling is the ratio of the correlation time tq to 
the thermal time-scale tt. At low temperatures, r < 1 (jr ^> tq), thermodynamics probes 
the FL exterior of the bubbles, but when r > 1 and r T <C tq, it reflects the QC interior of the 
bubbles. This accounts for the cross-over between FL and QC behaviors at r ~ 1. Moreover, 
T / B scaling over a wide range r ~ 1CT 1 to ~ 10 3 indicates that the quantum fluctuations in the 
ground-state are self-similar down to l/1000fh of the correlation time tq. 

The observation of zero-field quantum criticality in valence fluctuating /3-YbAlB 4 cannot 
naturally be interpreted as a conventional QCP, which would require a fortuitous combination of 
structure and chemistry to fine-tune the critical field B c to within 0.2 mT of zero. A more natural 
interpretation of the results is that /3-YbAlB 4 forms a quantum critical phase that is driven into a 
FL state by an infinitesimal magnetic field. The TIB scaling requires that the critical modes are 
Zeeman-split by a field, and as such, various scenarios, such as critical Fermi surfaces (28) or 
local quantum criticality with E/T scaling (6) may be possible contenders for the explanation, 
provided they can be stabilized as a phase. Established theoretical examples of a critical phase 
with T/B scaling include the Tomonaga-Luttinger liquid in half-integer spin chains and the 
one-dimensional Heisenberg ferromagnet (23, 24) . Experimentally, the d-electron metal MnSi 
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is a candidate for a quantum critical phase with anomalous transport exponents observed over a 
range of applied pressure (29). While present work provides a strong indication for existence of 
such a phase in /3-YbAlB 4 , future studies, in particular under pressure, are necessary in order 
to establish it definitively. 
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Figure 1 : (A) Crystal structures of (3- and a- YbAlB 4 , which are formed from straight and zigzag 
arrangements of distorted hexagons of Yb atoms (shaded in red for the a phase), respectively 
(IT). The crystallographic unit cells of both phases are orthorhombic and can be viewed as an 
interleaving of planar B-nets and Yb/Al-layers. (B) Magnetic part (/-electron contribution) of 
the specific heat Cm plotted as Cm/T versus T for both (3- (solid circles) and «-YbAlB 4 (open 
squares) (12). Cm/T at B = for the (3 phase shows a InT dependence for 0.2 K < T < 20 
K. T ~ 200 K was determined from the fit to C M /T = S /T In (T /T). The upturn in the 
lowest T may contain a nuclear contribution. Inset: Cm/T scaled by T compared with quantum 
critical systems CeCu 5 . 9 Au .i (T = 6.2 K) (14) and YbRhatSii^Ge^ (To = 24 K) (2, 16). The 
InT dependence of the three quantum critical materials collapse on top of each other using 
nearly the same coefficient So ~ 4 J/mole K, indicating a common meaning of T as the T scale 
below which ~ 70 % of the ground doublet entropy, R In 2, is released. 
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Figure 2: (A) Temperature dependence of the magnetic susceptibility M/B of both f3- (solid 
circles) and «-YbAlB 4 (open squares). The Curie- Weiss fit above 150 K yields a Weiss tem- 
perature 9w ~ — 110 K and an effective moment of ~ 2.2/ig for both systems. Inset shows the 
quantum critical B-T range where the scaling applies (solid circles in the blue shaded region) 
and the superconducting (SC) phase under the upper critical fields (open circles and triangles, 
see (12)). (B) Scaling observed for the magnetization at T IS 3 K and B 5 2 T. The data was fit- 
ted to the empirical Eq. 1 with scaling function <p{x) = Ax(A + x 2 )~ n , a form chosen to satisfy 
the appropriate limiting behavior in the Fermi liquid regime (72). The right inset shows Pear- 
son's correlation coefficient R for the fit with finite B c . Note that R reaches a maximum value of 
1 if the fit quality is perfect. The best fit is obtained with n — 1.25 ± 0.01 and B c = —0.1 ±0.1 
mT (light blue line), corresponding to a = 3/2 in the scaling form of the free energy, Eq. 2, 
see (12), and \B C \ < 0.2 mT. The left inset shows the B-T phase diagram of /3-YbAlB 4 in the 
low T and B region. The filled circles are determined from the peak temperatures of —dM/dT 
below which the FL ground state is stabilized. At low field, the thermodynamic boundary be- 
tween the FL and NFL regions is on a ksT ~ g^B line (broken line). The open circles are 
the temperature scale T FL below which the T 2 dependence of the resistivity is observed (9). 
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Supporting Online Material 

Materials and methods 

High purity single crystals of /3-YbAlB 4 were grown by a flux method (SI). Energy dispersive 
x-ray analysis found no impurity phases, no inhomogeneities and a ratio Yb:Al of 1:1. Surface 
impurities were carefully removed with dilute nitric acid before the measurements. A large 
mean free path exceeding 1000 A obtained from Shubnikov-de Haas measurements strongly 
supports the high purity of the samples (52). 

The magnetization data at T < 4 K and B < 0.05 T were obtained by using a high precision 
SQUID magnetometer installed in a 3 He- 4 He dilution refrigerator (S3). For fields above 0.05 T, 
we utilized a high precision Faraday magnetometer which is also installed in a 3 He- 4 He dilution 
refrigerator (S4). These magnetometers have a resolution of ~ 10~ 8 emu for the SQUID and 
~ 10~ 5 emu for the Faraday system. The samples used for the SQUID magnetometer measure- 
ments were high-purity single crystals with the residual resistivity ratio (RRR) > 200 (~ 30 
pieces, 0.82 mg). These samples were mounted using silver paste and were inserted in a pick- 
up coil and cooled by a heat link made of silver foils. The superconducting magnet was covered 
with a Nb superconducting shield and a ^-metal tube in order to eliminate the Earth's magnetic 
field. The pick-up and primary lines covered with Pb superconducting tubes were connected to 
a multiple purpose dc-SQUID probe located in a bath of liquid 4 He at 4.2 K. The magnetization 
M and ac-susceptibility Xac were obtained as dc and ac signals of the SQUID output. Xac below 
1 mT did not show any field dependence in the temperature range T > 0.08 K. This allows 
us to use the Xac data obtained at 0.31 mT and 0.62 mT to estimate the magnetization at the 
corresponding fields in the scaling analysis. For the Xac measurement, an ac field of 0.1 jll and 
frequency / = 16 Hz is applied along the c-axis. The residual magnetic field was estimated to 
be ~ 1.1 fiT using the dc-diamagnetic signals of the superconductivity of /3-YbAlB 4 measured 
under various magnetic fields of the order of f/T. For the Faraday magnetometer measurements, 
we used single crystals of 7.5 mg, whose typical RRR is as high as 140. The absolute values 
of M for both measurements were calibrated by comparing the data with those measured by a 
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commercial SQUID magnetometer at T > 2.0 K. 

The specific heat C was measured by a relaxation method. For the measurements above 0.4 
K, high-purity single crystals of 0.8 mg with RRR > 200 were measured by using a physical 
property measurement system (PPMS). The data below 0.4 K were obtained for high-purity 
single crystals of 2.2 mg with RRR > 200 by using a heat capacity cell installed in a 3 He- 4 He 
dilution refrigerator. The magnetic part of the specific heat Cm was obtained by subtracting the 
specific heat of a-LuAlB 4 , which is the non-magnetic isostructural counterpart of a-YbAlB 4 
(SI). To estimate C at B < 0.07 T in the analysis of Y H , a linear interpolation of the specific 
heat data at B — and 0.07 T was used (Fig. IB in the main text). The improved crystal 
quality, measurement technique, and more careful procedure to remove surface impurities of 
the crystals allowed us to measure more precisely the values of the low-temperature specific 
heat coefficient (e.g. C/T ~ 130 mJ/mol-K 2 at T = 0.4 K), which are fully reproducible 
and should replace the somewhat inaccurate values reported earlier in Ref. S5. Actually, in 
thermodynamic measurements in Ref. S5, a batch of about 50 crystals had to be used to gain 
signal intensity, with typical RRR values between ~10 and 300. 

Supporting online text 

1. Temperature dependence of the magnetizaiton and magnetocaloric ratio 

In order to characterize the power law behavior in the temperature dependence of the magneti- 
zation, we show —dM/dT versus T for selected fields along the c-axis on a logarithmic scale in 
Figure SI for both (3- (solid circle) and «-YbAlB 4 (open square). Below T ~ 3 K, —dM/dT for 
the (3 phase shows T^ 1 5 dependence, indicating the divergence of the susceptibility x ~ T -1 / 2 
(broken line) in the non-Fermi liquid region at T > B, and T-linear behavior (solid line) as 
expected for a Fermi liquid at T < B. The derivative of the magnetization —dM/dT under 
B = 6 T for the (3 -phase and under B — 0.1 T for the a-phase both show the peak at T* ~ 8 
K, due to the crossover from the high temperature local moment behavior to the Fermi liquid 
behavior at low temperatures. 
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The magnetocaloric ratio is defined as 



1 dS/dB 
TdS/dT 



OM/dT 
C H 



(S-l) 



The divergence of this quantity as T — > is a strong indicator of quantum criticality (S6). Our 
results for /3-YbAlB 4 in Fig. 3 in the main text show a clear divergence of Th/B as T — > 
in the non-Fermi liquid regime (B < T), whereas it levels off at low T in the Fermi liquid 
regime (T < B). Note that a small nuclear Schottky contribution to C (see Fig. IB in the main 
text) may slightly affect the results below ~ 0.2 K. However, after a subtraction of this nuclear 
contribution, the T dependence of Th becomes more divergent at low fields and temperatures 
than the results shown in Fig. 3 in the main text. In addition, according to theory, Th at a 
finite temperature is expected to exhibit two types of behavior as a function of field: in the NFL 
regime (B < T), T H oc \B - B c \ and in the FL regime (B > T), then T H oc 1/\B - B c \ (S6). 
Figure S2 shows the field dependence of T H , in which we see two clear regimes, as predicted 
by theory (S6). From the NFL regime, we can extract the critical field \B C \ < 0.2 mT, consistent 
with the estimate of B c obtained from the scaling behavior of M, Eq. 1 in the main text. 

2. Power-law fit to the low-temperature susceptibility 

The previous paper of /3-YbAlB 4 , Ref. S5, has initially reported the T -1 / 3 power-law of the 
low-temperature magnetic susceptibility. However present measurements, carried out on sig- 
nificantly higher quality single crystals, have shown that the exponent T" 1 / 2 provides a better fit 
to the experimental results and is also in agreement with the exponent derived from T/B scal- 
ing of magnetization (see next section). The present data more accurately represent the intrinsic 
behavior in the material, not only because of the higher quality of the single crystals used, but 
also because greater care was taken in the etching process to fully remove sample surface impu- 
rities, and furthermore, the current data were obtained without a drift in the SQUID output and 
without a sizable background signal. These refinements have been already reported in Ref. S3. 



18 



3. (T/B) scaling of thermodynamic quantities 

The scaling property of magnetization M shown in Eq. 1 of the main text, 

- dMjdT = B a ~ 2 (f)(T/B), (S-2) 

with a = 3/2 indicates that all thermodynamic properties can be expressed as a function of the 
ratio of T and B. Integrating both parts over temperature results in the following expression for 
magnetization: 

M c = M — M = —B a ~ 1 (f)(T/B), (S-3) 

where d(j)/dx = 4>{x). Above, M c and M is the critical and non-critical components of the 
magnetization, respectively. For a = 3/2, this means that M c / B 1 / 2 is a universal function of 
the ratio of (T/B) only, as plotted in Fig. S3. Since the magnetization M is the derivative of 
the free energy, M = —dF/dB, integrating both parts of Eq. IS -3 1 over field, one obtains the 
following simple scale-invariant form of the low-temperature free energy of /3-YbAlB 4 : 

F = B a f ( ^ ) , (S-4) 



B, 

where f(x) is a scaling function of the ratio T/B, which is related to <\> in Eq. lS-2l as follows: 
4>(x) = (a — l)f'(x) — xf"(x). To ensure that the free energy depends only on temperature in 
the non-Fermi liquid limit T ^> B, we require that the function f(x) ~ x a for x ^> 1. In the 
opposite limit T <^ B the system is a Fermi liquid, and hence one must be able to expand the 
free energy in powers of T 2 , which in turn requires that f(x) ~ const + 0(x 2 ) for x 1. In 
other words, f(x) is required to have the following asymptotic behavior in the two limits: 

f( . J x Q , forx3>l (T 3> B, non-Fermi liquid) 

* W ^ \ const + x 2 , for x < 1 (T < B, Fermi liquid) ( } 

The function chosen for the study is f(x) = — \(A + x 2 ) a ^ 2 , resulting in the scaling function 

<f)(x) = XAa(2 - a)x(A + x 2 )^' 2 = Ax(A + x 2 )~ n . (S-6) 

Figure 2B in the main text shows the fit of the empirical scaling relation Eq. IS-21 to the above 
form. The best fit is achieved with n = 1.25 ± 0.01 corresponding to a = 1.50 ± 0.02, in 
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agreement with the experimental data. This now allows us to write down the free energy in the 
simple form 

F = -^^{^Bf + ^TfY'\ (S-7) 

where g is the effective g-factor of the quantum critical excitations, related to the coefficient A in 
Eq. |S-6l as A = (g yU B / 'k B ) 2 , so that the scaling fit results in the effective moment g/j, s pa 1.94// B . 
The constant energy scale T in the prefactor can be determined from the following identity: 
T = (k B Aa(2—a)/A)'^ : , yielding k B T pa 6.56eV. While seemingly very large (of the order of 
the conduction electron bandwidth), this value of T agrees well with the measured value of the 
specific heat coefficient. Indeed, differentiating Eq. lS-7l twice over T, one arrives at the specific 
heat coefficient C QC /T ~ ^p(TT)" 1 / 2 , which at T = 0.1 K yields a value 77 mJ/mol K 2 for 
the quantum critical component of the specific heat coefficient. This value of Cqq/T is close to 
an experimental estimate of the QC component of C/T at T — 0.1 K, 7 QC ~ 60 mJ/mol K 2 for 
the (3 phase. The estimate is obtained after subtracting C/T pa 130 mJ/mol K 2 for ct-YbAlB 4 
at 0.1 K, the sum of the electronic non-critical component and the nuclear Schottky component, 
both of which are likely similar for these locally isostructural compounds. 

The above simple form of the scaling relation (IS-71) is significant in that the quantum critical 
free energy only depends on the distance from the origin in the (k B T, g^ B B) coordinates, similar 
to e.g. the case of the Tomonaga-Luttinger liquid, where the free energy takes a similar form 
F TL oc (k B T) 2 + (gfi B B) 2 with a model-dependent ^-factor g (S10). 

The above scaling relations, Eqs. lS-2[|S-4l can be easily modified to a more general case of 
non-zero critical field by setting B — > B — B c in the above expression. This allowed us to place 
an upper bound on the value of the critical field, \B C \ < 0.2 mT, deduced from Pearson's fitting 
quality shown in the inset of Fig. 2B, as discussed in the main text. 

The field dependence of the specific heat coefficient j(T) = C(T) /T = —d 2 F/dT 2 can be 
deduced from Eq. lS-4l to be 

7(T,£) = -B a ' 2 f"(T/B). (S-8) 
The limiting T — > behavior of this quantity is proportional to the effective mass of the quasi- 
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particles: m*(B) oc j(T = 0, B), giving 



m*{B) oc B a - 2 f"{0) 



(S-9) 



which yields m* oc B^ 1 ^ 2 for the empirically determined value of a = 1.5. 

T/£> scaling in /3-YbAlB 4 is observed over an unprecedentedly wide range of more than 
three decades in the non-Fermi liquid regime T/B > 1. This may be contrasted with other 
known examples of T/\B — B c \ scaling in quantum critical heavy fermion materials such as 
CeCu 6 _ x Au x (B c = 0) and YbRh 2 Si 2 (B c ^ 0), where in the non-Fermi liquid regime, the 
scaling was observed over roughly one decade of this ratio (S7,S9). Known theoretical examples 
of quantum models displaying T/B scaling include (S10): the Tomonaga-Luttinger liquid in 
one-dimensional metals, the critical (2+l)-dimensional 0(N) sigma model and a drained Fermi 
liquid with a chemical potential fi = 0. 

Currently, it is not possible to extract the pure electronic contribution to the specific heat at 
zero field because the nuclear contribution is unknown. However, we may analyze the change 
in the specific heat under application of a low field using the Maxwell relation (dS/dB)x = 
(dM/dT) B and the scaling equation (IS-21) . Based on the Maxwell relation, the specific heat 
C(B) in magnetic field B can be expressed as follows, 



Here, C(0) is the zero field specific heat. Figure S4 shows that the calculated results using 

Eq. IS-21 and Eq. IS- 1 01 are in the full agreement with the experimental results for B = 0.07 
and 0.5 T. This indicates that the free energy given in Eq. 2 in the main text and Eq. IS-41 is 
consistent with the temperature and field dependence of the specific heat in the quantum critical 
regime. 




(S-10) 
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1 0" 1 10° 1 o 1 

r (K) 

Figure SI: —dM/dT versus T for selected fields on a logarithmic scale for both 0- (solid 
circles) and «-YbAlB 4 (open squares). The sudden downturn below 0.1 K in the low field data 
corresponds to the onset of superconductivity. 
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B (T) 

Figure S2: Magnetic field dependence of the magnetocaloric effect for /3-YbAlB 4 , T H = 
-(dM/dT)/C, obtained from -dM/dT and the total specific heat C. In the NFL regime 
at T > B, T H increases linearly with B, while in the FL regime at T < B, it exhibits B~ l de- 
pendence, consistent with theory (S6). The upturn observed at ~ 2 T in the temperature range 
of 0.4 K < T < 1 .5 K indicates the boundary between the quantum critical and non-quantum 
critical regions. 
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Figure S3: Scaling plot of Mjy/B vs. ratio T/B for /3-YbAlB 4 , which according to Eq. dO 
should be a universal function of this ratio. Here, M c is the critical contribution to the magne- 
tization: M c = M — M , where we have subtracted M = XoB using the constant non-critical 
susceptibility xo — 0.017 emu/mol. This latter constant term, whose nature is immaterial to 
the scaling analysis, is close to the zero-T susceptibility of the non-critical a-YbAlB 4 . It may 
originate from constant Van Vleck contribution to susceptibility and/or from Pauli susceptibil- 
ity of the non-critical parts of the Fermi surfaces. Such a subtraction is well justified in the 
literature and was used for instance to establish the T/B scaling in CqC\Xq_ x Pm x (see Fig. 4b 
in Ref. (ST)). Under fields B > 0.5 T, M c does not follow a single scaling curve, because 
of a small field-nonlinear contribution to M , which is neglected in the analysis. Inset shows 
the quantum critical B-T range where the scaling applies (solid circles in the blue shaded re- 
gion) and the superconducting phase boundary (open circles (S3) and open triangles (S8)) of 
/3-YbAlB 4 . 
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T(K) 

Figure S4: Specific heat divided by temperature C/T in the low fields below 0.5 T along 
the c axis (filled symbols) and those estimated by using the scaling relation, Eq. (1) in the 
main text (or Eq. dS-21) ) as well as the Maxwell relation, Eq. (|S-10I) (solid lines, see text). The 
experimental results and the calculation agree well with each other in the low field below 0.5 T 
in the quantum critical regime. The broken line is a fit to the zero filed data C(0) /T, which was 
used in Eq. (fS^TOl . 
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